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Abstract 

In [5] the authors obtained Mean-Field backward stochastic differential equations (BSDE) 
associated with a Mean-field stochastic differential equation (SDE) in a natural way as limit of 
r — some highly dimensional system of forward and backward SDEs, corresponding to a large number 

^vO • of "particles" (or "agents"). The objective of the present paper is to deepen the investigation of 

such Mean-Field BSDEs by studying them in a more general framework, with general driver, and 
to discuss comparison results for them. In a second step we are interested in partial differential 
equations (PDE) whose solutions can be stochastically interpreted in terms of Mean-Field BSDEs. 
For this we study a Mean-Field BSDE in a Markovian framework, associated with a Mean-Field 
forward equation. By combining classical BSDE methods, in particular that of "backward semi- 
groups" introduced by Peng [14] . with specific arguments for Mean-Field BSDEs we prove that 
this Mean-Field BSDE describes the viscosity solution of a nonlocal PDE. The uniqueness of this 
viscosity solution is obtained for the space of continuous functions with polynomial growth. With 
the help of an example it is shown that for the nonlocal PDEs associated to Mean-Field BSDEs 
one cannot expect to have uniqueness in a larger space of continuous functions. 
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1 Introduction 

Classical mathematical mean-field approaches can be met in different fields, so in Statistical 
Mechanics and Physics (for instance, the derivation of Boltzmann or Vlasov equations in the kinetic 
gas theory) and in Quantum Mechanics and Quantum Chemistry (for instance, the density func- 
tional models or also Hartree and Hartree-Fock type models). In a recent series of papers (see |10j 
and their papers cited therein) Lasry and Lions extended the field of such mean-field approaches 
to problems in Economics and Finance, and also to the theory of stochastic differential games. 
Recently Buckdahn, Li and Peng [5] studied a special mean field problem in a purely stochastic 
approach and deduced a new kind of BSDEs which they called Mean-Field BSDEs. In order to be 
more precise, we consider the stochastic dynamics X = (Xt)te[o,T\ of a particle 

dX t = a(t,X t ,X},...,Xl f )dB t + b(t,X t ,Xl...,X?)dt,t€[0,T\, 

X = x G R n , { ' ' 

which is governed by a d-dimensional Brownian motion B = (-£>*) te[o,T] and influenced by the 
positions X , . . . ,X of N other particles that are supposed to be mutually independent of each 
other, also independent of the driving Brownian motion B, and of the same law as X. The time 
T > is an arbitrarily fixed horizon. In [5] it is shown that, under appropriate assumptions, SDE 
(1.1) has a unique solution. 

For the study of the limit behavior N — ► +oo of the above SDE of rank N, the authors of [5] 
assumed that the coefficients a = a N and b = b N of SDE (1.1) of rank N have the form 

1 N 1 N 

a N (t,x,x\...,x N ) = — ^a(i,x\x), b N (t,x,x 1 , . . . ,x N ) = — J^ b(t, x l ,x), 

i=l i=l 

(t,x, (s 1 , . . .,x N )) G [0,T] x R n x R nxN . The functions a : [0,T] x R n x R n -> R nxd , b : [0,T] x 
R n x R n — > R n are supposed to be bounded, measurable and, moreover, Lipschitz in the space 
variables, uniformly with respect to the time variable. Denoting by (X( N > , J^' 1 , . . . ,X^ N '' N ) := 
(X, X 1 , . . . ,X N ) the solution of SDE (1.1) with the coefficients a N and b N , it is shown in [5j 
that the the processes l'") converge uniformly on the time interval [0,T], in L 2 , to the unique 
continuous F B -adapted solution X of the limit equation (Mean-Field SDE) 



dX f = E 



a(i,X t ,/x) . dB t + E a(t,X t ,n) . dt, t G [0,T], X = x (1.2) 



n=x t 

(¥ B is the filtration generated by B). 

With the SDE of rank N the authors of [5j associated a BSDE of the form 

dY t = -f N (t, X t ,Xl ..., X t N ,Y t , Y t \..., Y t N , Z t , Z\, . . . , Z t N )dt + Z t dB t , t G [0, T], 
Y T = $ N (X T ,Xi,...,X»), 



(1.3) 



in which, in the same spirit as for the forward SDE of rank N, (X*, Y l , Z l ), 1 < i < N, is a sample 
of N mutually independent triplets which are independent of the driving Browning motion B and 
obey the same law as the solution (X, Y, Z); recall that (X, X 1 , . . . , X ) is the solution of the above 
SDE (1.1) with coefficients a and b . The coefficients / , $ for the above BSDE are introduced 
in the same spirit as a , b N : 
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/ (t,x,x ,...,x ,y,y ,...,y ,z,z,...,z ) = —^2 f{t,x\x,y\y, z l , z), 

i=l 

1 N 

$ N (x,x 1 ,...,x N ) = -Y,Hx\x), 



N 

i=l 



te [0,T], x,xi £R n ,y,y j £R,z,z^ e R d , 1 < j < N, where / : [0,T]xR n xR n xRxRxR d xR d - 
R and <I> : i? n x i? n — > R are bounded, measurable and Lipschitz in the variables (x', x, y' , y, z', z) 
uniformly with respect to t € [0, T]. In [5] it is proved that, with 



N\ 



(f (at )j Y w,l t Y (N), Nj Z (N) )Z W,\ Z (N), N) :={Y ,Y\...,Y N ,Z,Z\...,Z 

denoting the solution of the above BSDE (1.3), the couple of processes (Y^ \Z^ '), interpreted 
as random variable with values in C([0,T];R) x L 2 ([0,T];R d ), converges to the unique square 
integrable, F s -progressively measurable solution of the limit equation (called Mean-Field BSDE) 



dY t = -E 
Y T = E 



f(t,X t ,n,Y u \,Zt,Q . . . dt + ZtdB u te [0,T], 

lf l =Xt,X=Y t ,(=Zt n^\ 



*(Xr,/x) 



ti=X 7 



where X is the solution of Mean-Field SDE (1.2). 

The objective of the present paper is to deepen the investigation of the above Mean-Field 
BSDE. In a first leg we study the existence and uniqueness for Mean-Field BSDEs in a rather general 
setting, with drivers which, in difference to [5], are not necessarily deterministic coefficients. In 
addition to the existence and the uniqueness also the comparison principle for this new type of 
BSDEs is discussed and illustrated by examples. 

The main objective of the paper concerns the study of Mean-Field problems in a Markovian 
setting. To be more precise, we investigate Mean-Field BSDEs associated with Mean-Field forward 
SDEs and partial differential equations (PDEs) whose solutions are described by them. The system 
dynamics we investigate is given by the following SDE 



j dXt't = E[b(s,X° s ' x ° ,/x)] ,, c ds + E[a(s,X^ x ° ,»)] t ,<dB a , s £ [t,T], 

I xl'< = c 

Precise assumptions on the coefficients b : [0, T] x R n x R n -*• E n and a : [0, T] x R n x R n 
are given in the following sections. 

With the above SDE we associate the BSDE: 

J -dY l s ' C = E[f(s, X^ , /x, Y S °' X0 ,\, C)] M=X « x=Y t,c (=zl xds - Z f /dB s , s € [t, T], 

{ y« = e^x ™^)}^. 



(1.5) 



pnxd 



(1.6) 



Under the assumptions on / and <£ that are introduced in Section 5, the above BSDE has a unique 
solution (Y s ' x , Z s ' x ) s <=it i T] an d we can define the deterministic function 



u(t,x) = Y t l > x . (1.7) 



We prove that u(t,x) is the unique viscosity solution in C p ([0,T] x R n ) to the following nonlocal 
PDE 

§- t u(t, x) + Au(t, x) + E[f(t, Xt' xo ,x, u(t, X°' xo ),u(t, x), Du(t, x).E[a(t, X^ xo , x)])] = 0, 

(t,x) G [0,T) xl", 



u(T,x) =E[$(Xt ,x °,x)}, xe 



(I- 



with 



Au(t, x) := -tr(E[a(t, X^ x ° , x)]E[a(t, X^ x ° , x)] T D 2 u(t, x)) + Du(t, x).E[b(t, X^ x ° , x)]. 

In particular, it is shown that the space C p ([0, T] x M. n ) is the optimal space in which the uniqueness 
can be got. 

Our paper is organized as follows. Section 2 recalls some elements of the theory of BSDEs 
which are needed in what follows. Section 3 investigates the properties of general Mean-Field 
BSDEs. We first prove the uniqueness and existence of the solution of Mean-Field BSDE (The- 
orem 3.1) but also the comparison theorem (Theorem 3.2) and the converse comparison theorem 
(Theorem 3.3). Similarly we investigate Mean-Field Forward SDEs in Section 4. In Section 5 we 
investigate decoupled Mean-Field Forward-Backward SDEs (FBSDEs). Their value function u (see 
(5.4)) turns out to be a deterministic function which is Lipschitz in x (see (5.5)) and ^-Holder con- 
tinuous in t (Theorem 5.2). Moreover, it satisfies the dynamic programming principle (DPP) (see 
5.10). A key element in the proof of the DPP is the use of Peng's backward semigroups (see (14J). 
We change slightly its definition; this allows to shorten the argument for the proof that u is a 
viscosity solution of the associated PDE (Theorem 6.1). Finally, the uniqueness of the viscosity 
solution in the space of continuous functions with polynomial growth is proved in Section 7. 

2 Preliminaries 

Let {Bt}t>o be a d-dimensional standard Brownian motion defined over some complete prob- 
ability space (£l,J-, P). By F = {J- s , < s < T} we denote the natural filtration generated by 
{-B s }o< s <T and augmented by all P-null sets, i.e., 

T s = a{B r ,r< s} VA/>, a G [0,T], 

where Mp is the set of all P-null subsets and T > a fixed real time horizon. For any n > 1, 
\z\ denotes the Euclidean norm of z G W 1 . We also shall introduce the following both spaces of 
processes which are used frequently in what follows: 

iSf(0,T;R) := {(V>t)o<t<T real-valued F-adapted cadlag process : 

E[ sup |V>t| 2 ] < +cx)}; 
0<t<T 

W|(0,T;]R n ) := {(ipt)o<t<T M n -valued F-progressively measurable process : 

IIV<ll!=MJo T |^| 2 ^]<+cx)}- 

Let us now consider a function g : Q x [0, T] x R x Mr — > R with the property that 
(g(t,y,z))t£[o t T\ i s F-progressively measurable for each (y,z) in R x R , and which is assumed 
to satisfy the following standard assumptions throughout the paper: 
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(Al) There exists a constant C > such that, dtdP-a.e., for all 2/1,2/2 £ R ; ^l^ £ R ; 
|c/(t,yi,zi) - g(t,y 2 ,z 2 )\ < C(\yi -y 2 \ + |zi -^|); 

(A2) 5 (-,0,0)eH|(0,T;R). 

The following result on BSDEs is by now well known, for its proof the reader is referred to 
Pardoux and Peng [TZ] . 

Lemma 2.1. Under the assumptions (Al) and (A2), for any random variable £ € L 2 {@,J-t, P), 
the BSDE 

Vt = £+ I g(s,y s ,z s )ds- I z s dB s , 0<t<T, (2.1) 

Jt Jt 

has a unique adapted solution 

(yuz t ) t£[0tT] eS$(0,T;R)xH 2 ¥ (0,T;R d ). 

In what follows, we always assume that the driving coefficient g of a BSDE satisfies (Al) and 
(A2). 

Let us remark that Lemma 2.1 remains true when assumption (Al) is replaced by weaker 
assumptions, for instance those studied in Bahlali [lj, Bahlali, Essaky, Hassani and Pardoux [2] or 
Pardoux and Peng [T3]. However, here, for the sake of simplicity of the calculation we prefer to 
work with the Lipschitz assumption. 

We also shall recall the following both basic results on BSDEs. We begin with the well-known 
comparison theorem (see Theorem 2.2 in El Karoui, Peng and Quenez [7]). 

Lemma 2.2. (Comparison Theorem) Given two coefficients g\ and g 2 satisfying (Al) and (A2) 
and two terminal values £1, £2 € L 2 (0, Ft, P), we denote by (y 1 ,^ 1 ) and (y 2 ,z 2 ) the solution of 
the BSDE with data (£1,51) and (£2,52); respectively. Then we have: 

(i) (Monotonicity) If £1 > £ 2 and g\ > g 2 , as., then y\ > y 2 , a.s., for all t £ [0,T]. 

(ii) (Strict Monotonicity) If, in addition to (i), we also assume that P(£i > £2) > 0, then 
P{lJt > Vt} > 0) <t <T, and in particular, y$ > y 2 . 

Using the notation introduced in Lemma 2.2 we now suppose that, for some g : f2 x [0, T] x 
IxK 1 * — ► R satisfying (Al) and (A2), the drivers gi, i = 1, 2, are of the form 

9i(s,yl,zl) = g(s,y\,z\) +<Pi(s), dsdP-a.e., 

where (fi S 7i|(0,T;R). Then, for terminal values £1, £2 belonging to L 2 (S7, Tt^P) we have the 
following 

Lemma 2.3. The difference of the solutions {y l ,z l ) and (y 2 ,z 2 ) of the BSDE with data (£1,51) 
and (£,2,92), respectively, satisfies the following estimate: 

la 1 - v 2 t? + kE\gS^U - vl\ 2 + 1*1 - *2I>I*] 

< S[^ T -«|6-6| 2 |^t]+^LA r ^«-*)|y, 1 ( g )-y, 2 ( 8 )| 2 ds|^] J P-a.s., forallO<t<T, 

where (3= 16(1 + C 2 ). 

For the proof the reader is referred to Proposition 2.1 in El Karoui, Peng and Quenez [7] or 
Theorem 2.3 in Peng [14]. 
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3 Mean-Field Backward Stochastic Differential Equations 

This section is devoted to the study of a new type of BSDEs, the so called Mean-Field BSDEs. 

Let (Cl,F,P) = (fixfl,f®f,P®P) be the (non-completed) product of (n,F,P) with 
itself. We endow this product space with the filtration F = \Tt = T <£> J-t, < £ < T}. A random 
variable £ G L (Q, T , P;W l ) originally defined on Q is extended canonically to O: £'(u/,u;) = 
£(u/), (u',u) G Cl = n x Q. For any 9 G L 1 ((i,f',P) the variable 0(.,u>) : -> R belongs to 
L l {Q,T,P), P(duj)-a.s.; we denote its expectation by 



£'[0(.,w)] = / 6{tJ ,u)P{dJ). 
Jci 

Notice that £'[#] = E'[9{.,oj)} G L 1 ^,?^), and 

£[#](= / fidP = / £'[0(.,u;)]P(du;)) = £[£'[#]]• 

The driver of our mean-field BSDE is a function / = f(u',to,t,y',z f ,y,z) : Cl x [0,T] x R x R d x 
R x R — > R which is F-progressively measurable, for all (y',z',y,z), and satisfies the following 
assumptions: 

(A3) There exists a constant C > such that, P-a.s., for all t € [0, T], yi,y 2 ,y[,y 2 € 

R, 31,32,4,4 GR rf , 

l/(*,yi,4>yi>2i) - f(t,y^ z ^ 2/2,^2)1 < c(bi -y 2 l + 14 -41 + I2/1 — 2/2I + ki -^D- 

(A4) /(-,(), 0,0,0) GH?(0,T;R). 

Remark 3.1. Ze£ (5 : £1 x [0, T] — > R, 7 : f2 x [0, T] — > M. d be two square integrable, jointly 
measurable processes. Then, for our driver, we can define, for all (y,z) G R x R rf , dtP(dto)-a.e., 

fP>t(u,t,y,z) = E'[f(.,iv,t,P' tn ' t ,y,z)} 

= f n f(Lu',co,t,/3t(u>'),7t(u'),y,z)P(dw'). 

Indeed, we remark that, for all (y,z), due to our assumptions on the driver f, (f(.,t,(3' t ,'j' t ,y,z)) G 
Wl(0, T;R), and thus f /3 '' y (.,.,y,z) G 7Y|(0,T;R). Moreover, with the constant C of assumption 
(A3), for all ( Vl , Zl ), (y 2 , z 2 ) G R x R d , dtP(du)-a.e., 

\f^(u,t, yi , Zl ) - f^(u>,t,y 2 ,z 2 )\ < C(\ yi -y 2 \ + \ Zl - z 2 \). 

Consequently, there is an ¥ -progressively measurable version of /^' 7 (., -,y, z), (y,z) G R x M. d , 
such that /^' 7 (o;, t, ., .) is dtP{duo)-a.e. defined and Lipschitz in (y, z); its Lipschitz constant is that 
introduced in (A3). 

We now can state the main result of this section. 

Theorem 3.1. Under the assumptions (A3) and (A4), for any random variable £ G L 2 (0, JF r , P), 
the mean-field BSDE 

Yt = £+ f E'lf(sX,Z' s ,Y s ,Z s )]ds- f Z s dB s , 0<t<T, (3.1) 

Jt Jt 

has a unique adapted solution 

(Y t ,Z t ) te[0jT] eS§(0,T;R)xTiU0,T;R d ). 
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Remark 3.2. We emphasize that, due to our notations, the driving coefficient of (3.1) has to be 
interpreted as follows 

E'[f(sX,Z' s ,Y s ,Z s )}(co) = E'[f(8,Y!,Zs,Y s (u),Z s (w)] 

= J a f(u/,w,8,Y a (cS),Z s (u/),Y s (u;),Z s (u>))P(du/). 



Proof. We first introduce a norm on the space 7Yp(0, T; R x R ) which is equivalent to the canonical 

norm: 

T ,h 



\H-)\b = \Ej KlV'dsj , /3>0. 

The parameter /3 will be specified later. 

Step 1 : For any (y,z) G W|(0,T;R x R d ) there exists a unique solution (Y,Z) G S^(0,T;\ 

HJ(0,T;R d ) to the following BSDE: 

Y t = Z+ [ E'[f(s,y' s ,z' s ,Y s ,Z s )]ds- f Z s dB s , 0<t<T. 
Jt Jt 



(3.2) 



Indeed, we define g( y,z >(s, [i, v) = E'[f(s, y' s , z' s , fi, v)\. Then, due to Remark 3.1, g( y,z >(s, /i, v) 
satisfies (Al) and (A2), and from Lemma 2.1 we know there exists a unique solution (Y, Z) G 
<S 2 (0, T;R) x H 2 (0,T;R d ) to the BSDE (3.2). 

Step 2: The result of Step 1 allows to introduce the mapping (Y.,Z.) = I[(y.,z.)] : 7f|(0,T;R x 
Rrf) -> Hj(0, T;lx R d ) by the equation 

Y t = i+ I E'[f(s,y' s ,z' s ,Y s ,Z s )]ds- f Z s dB s , < t < T. (3.3) 

Jt Jt 



For any (y\z v ), (y 2 ,z 2 ) G H 2 (0,T;R x R d ) we put (Y 1 ^ 1 ) = I[{y\z v % (Y 2 ,Z 2 ) = I[(y 2 ,z 2 )}, 
(y,z) = (y 1 — y 2 ,z x — z 2 ) and (Y,Z) = (Y 1 - Y 2 ,Z X — Z 2 ). Then, by applying Ito's formula to 
e^ s \Y s \ 2 and by using that Y 1 , Y 2 G <S|(0, T; R) we get 

= E[tf e^-t) 2 Y r {g^^){r,Y r \Zl)-g^^){r,Y 2 ,Z 2 ))dr\T t i te [0,T]. 
From assumption (A3) we obtain 

(f - 2C- 2C 2 )S[r T e ^|i;| 2 (ir] + ^[J^e^l^r] 
< ^{S[/ T e^|ifr| 2 dr] + S[/ r e^|^.| 2 dr]}. 

Thus, taking /3 = 16C 2 + 4C + 1 we get 

£[/" e^(|Y r | 2 + |Z r | 2 )dr]<;W e^(|y r | 2 + |£ r | 2 )(ir], 
Jo 2 Jo 

that is, ||(Y, Z)^ < -^ 1 1 (y, 2)11/3- Consequently, I is a contraction on 7^.(0, T;R x ^ d ) endowed 
with the norm ||.|U, and from the contraction mapping theorem we know that there is a unique 
fixed point (Y, Z) G 7^(0, T; R x R d ) such that /(Y, Z) = (Y, Z). On the other hand, from Step 1 
we already know that if I(Y, Z) = (Y, Z) then (Y, Z) G <S|(0, T; R) x W|(0, T; R d ). D 
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Using the notation introduced in Theorem 3.1 we now suppose that, for some /:Ox[0,T]x 
IxR^xtxl^R satisfying (A3) and (A4), the drivers fc, i = 1, 2, are of the form 

f(s, Yf,Z^Y:, Z\) = f(s, Yf, Zi',Yi, Z\) + Vl (s), dsdP-a.e., i = 1, 2, 

where ifi G Hj{0, T; R) and (Y\ Z { ) is the solution of Mean-Field BSDE (3.1) with data (/*, &), i = 
1,2. Then, for arbitrary terminal values £i, £2 belonging to L 2 {Q,,J : t,P) we have the following 

Lemma 3.1. T/ie difference of the solutions (Y , Z ) and (Y 2 ,Z 2 ) of BSDE (3.1) with the data 
(£ii/i) anc ^ (£21/2), respectively, satisfies the following estimate: 

e[\y? - y??] + ^[jf^Hin 1 - n 2 i 2 + izi - ^i 2 )&] 

< ^[e^ T -')[ei -6I 2 ] +£[jfe^-*)|v9i(s) -^ 2 (s)| 2 ds], /or allO<t<T, 
where (3 = 16(1 + C 2 ). 

The proof uses a similar argument as that of the proof of Theorem 3.1 and is therefore 
omitted. 

Now we discuss the comparison principle for Mean-Field BSDE. We first give two counterex- 
amples to show that if the driver / depends on z' or / is decreasing with respect to y' we can't get 
the comparison theorem. 

Example 3.1 For d = 1 we consider the Mean-Field BSDE (3.1) with time horizon T = 1, with 
driver f(u>',uj,s,y',z',y,z) = —z' and two different terminal values £1, £2 S -^ 2 (^,-^T ; -P)- Let us 
denote the associated solutions by (Y 1 ^ 1 ) and (Y 2 ,Z 2 ), respectively. Then, 

Yt = & + / E[-Z\\ds - I ZldB s , < t < 1, i = 1, 2. (3.4) 

Jt Jt 

We let £1 = -(-B+) 3 and define Y? = Y/ + jf E[Z\]ds. Then (Y 1 ,^ 1 ) is the unique solution 
of the BSDE Y^ 1 = 6 - ^ Z\dB s , t G [0,T]. Thus, Yq 1 = £?[&] = -£[(£+) 3 ] = --|=. On 



the other hand, from the Clark-Ocone formula we know that Z 1 is the predictable projection of 
the Malliavin derivative (DtOte[o,T] of £ (Dt£ denotes the Malliavin derivative of £ at time i; the 
interested reader is referred, e.g., to Nualart [H]). This implies E\Z\\ = E[D t £\ = E[— 3(Bf) 2 ] = 
-§, i G [0, 1]. Therefore, Yq 1 = Y^ - ff E[Z}}ds = --^- + | > 0. Let now £ 2 = 0. Then, obviously, 
(Y 2 , Z 2 ) = (0,0). Hence, we have Yq 1 > Y 2 although £1 < &, P-a.s. and P{£i < 6} > 0. 
Example 3.2 Let again d = 1. We consider Mean-Field BSDE (3.1) driven by the function 
f(ui',uj,s,y',z',y,z) = —y', with time horizon T = 2 and two different terminal values £1, £2 £ 
L 2 (fL,T T ,P): 

17 = £ + f E[-Y l s ]ds - f Z'idB s , < t < 2, i = 1, 2. (3.5) 

Jt Jt 

By choosing first £1 = (i?i) 2 we have -EfYj 1 ] = e~^ 2_t \ t G [0,2]. Furthermore, for t G [1,2] we 
have Yl = (B1) 2 - ff e'^'^ds = (B1) 2 - (1 - e^ 2 "*)) and Z^ 1 = 0, P-a.s. Consequently, Y/ = 
(Bi) 2 — (1 — e _1 ) < on the set {(Si) 2 < 1 — e -1 } which is of strictly positive probability. Finally, 
for £2 = we have the solution (Y 2 ,Z 2 ) = (0,0). Therefore, in our example P(Y^ < Y 2 ) > 0, 
although £1 > £,2, P-a.s. 

The above both examples show that we cannot hope to have the comparison principle between 
two Mean-Field BSDEs whose drivers depend both on z' or are both decreasing in y' '. 



Theorem 3.2. (Comparison Theorem) Let fi = fi(u>,t,y' ,z' ,y,z), i = 1,2, be two drivers satisfy- 
ing the standard assumptions (A3) and (A4). Moreover, we suppose 
(i) One of both coefficients is independent of z' ; 
(ii) One of both coefficients is nondecreasing in y' . 

Let £1, £2 & L 2 (Q,J-t,P) and denote by (Y^Z 1 ) and (Y 2 ,Z 2 ) the solution of the mean-field 
BSDE (3.1) with data (£1, f\) and (£, 2 , f^), respectively. Then if^\ < £ 2 , P-a.s., and f\ < f'2, P-a.s., 
it holds that also Y/ < Y t 2 , t G [0,T], P-a.s. 

Remark 3.3. The conditions (i) and (ii) of Theorem 3.2 are, in particular, satisfied, if they hold 
for the same driver fj but also if (i) is satisfied by one driver and (ii) by the other one. 

Proof. Without loss of generality, we assume that (i) is satisfied by f\ and (ii) by f 2 . Then, with 
the notation (Y, Z) := (Y 1 - Y 2 , Z 1 - Z 2 ), £ := & - £ 2 , 

r T r T - 

Y t = i+ E'[h{s,Y}',Y},Zl)-f 2 {s,Y 2 ',Z 2 J,Y 2 ,Z 2 s )]ds- Z s dB s , 0<t<T, (3.6) 
Jt Jt 

and from Poo's formula applied to an appropriate approximation of f{y) = {y + ) 2 , y € R, in which 
we take after the limit, we obtain 

E[(Y t + ) 2 ] + E[f t T I { y s>0} \Z s \ 2 ds] 
= 2E[j t T Y+{E'[h{s, Y}',Y},Z\) - f 2 (s, Y 2 ',Z 2 ',Y 2 , Z 2 )])ds] 
< 2E[f t T Y+(E'[h{s, Y}'X,Z]) - Ms, Y}',Y 2 , Z 2 )])ds] 

+2E[f t T Y+(E'[f 2 (s, Y}',Z 2 ',Y 2 , Z 2 ) - f 2 {s, Y 2 ' , Z 2 ',Y 2 , Z 2 )))ds), t G [0, T]. 

Since f 2 (s,y' ,z',y, z) is nondecreasing in y' we get from (A3) 

E[Y+{E'[h{ S ,Y}' \Y},Z\) - h{s,Y}' \Y 2 ,Z 2 )])] 
+E[Y+(E'[f 2 (s, Y}\ Z 2 ',Y 2 , Z 2 ) - f 2 (s, Y 2 ' , Z 2 ',Y 2 , Z 2 )})} 

< E[Y+(E'[C(Y}' - Y 2 ') + + C\Y S \ + C\Z S \})\ 
= C(E[Y+\) 2 + CE[{Y+) 2 ] + CE[Y+\Z s \} 

< (2C + C 2 )E[(Y+) 2 ] + \E[I { y s>0} \Z s \ 2 ], se[0,T}. 

Consequently, 

E[(Y+) 2 \ + l -E[j^ I {Ys>0} \Z s \ 2 ds]<(4C + 2C 2 )J t E[(Y s + ) 2 ]ds, t E [0,T], 



from where we can conclude with the help of Gronwall's Lemma that Y^ — Y 2 = Yt < 0, t € 
[0,T], P-a.s. D 

We also have a converse comparison theorem. 

Theorem 3.3. (Converse Comparison Theorem) We retake the assumptions of Theorem 3.2 and 

suppose that, additionally, for some t € [0, T], Y^ = Y 2 , P-a.s. Then 

(i) Y} = Y 2 , s € [t,T], P-a.s., and 

(ii) If f 2 satisfies (ii) of Theorem 3.2 then B '[h(s ,Y}' , Z}' ' ,Y 2 , Z 2 )} = E'[f 2 (s,Y 2 ' , Z 2 ' ,Y 2 , Z 2 )}, 

dsdP-a.e. on [t,T], and if fi satisfies (ii) then the symmetric result holds. 



Proof. We use the notation introduced in the proof of Theorem 3.2 and suppose again that /i 
satisfies (i) and f 2 (ii). Then, from the Lipschitz property of j\, there exist some F-progressively 
measurable, bounded processes a, (3, defined over £1 x [0, T], such that 

h(s, Y}\y}, Z\) - f^s, Y}',Y 2 , Z 2 ) = a s Y s + (i s Z s , s G [0, T], 

and since Y 1 , Z 1 , Y 2 , Z 2 don't depend on a/, 

E'[h{s, Y}',Y}, Z]) - h{s, Y}',Y 2 , Z 2 )} = E'[a s ]Y s + E'%]Z S , s G [0, T]. 

Thus, from Ito's formula, 

Y t = exp{/ t T E>[a r ]dr}Z + jf exp{// E , [a r ]dr}E'[f 1 ( a , Y}' ,Y 2 , Z 2 ) - Ms, Y 2 ' , Z 2 ',Y 2 , Z 2 )]ds 
— J. exp{J t E'[a r ]dr}Z s dB s , P-a.s., 

where B s = B s — L E'[(3 r }dr, s G [0, T]. It is well known that B = (B s ) is an (F, P)-Brownian 
motion with P := exp{ f Q E'[(3 s \dB s — \j \E'[(3 S \ \ 2 ds}P. From the boundedness of a and (5 we then 
deduce easily that /. exp{J t E'[a r ]dr}Z s dB s is an (F, P)-martingale increment. Consequently, 

= Y t = E p [exp{f t u E>[a r ]dr}Y u \F t ] 

+E p [J t u exp{// E'[a r \dr}E'[his, Y^',Y 2 , Z 2 ) - / 2 ( S , Y 2 ' , Z 2 ',Y 2 , Z 2 )]ds\T t ], 

P-a.s., for all t < u < T. 

To conclude it suffices now to recall that, due to Theorem 3.2, Y u = Y^ — Y 2 < 0, P-a.s. and 

f\is,Y}',Y 2 ,Z 2 s ) < f 2 is,Y}',Z 2 ',Y 2 ,Z 2 ) < his,Y 2 ',Z 2 ',Y 2 ,Z 2 ), dsdP-a.e. on Q x [t,T]. 

D 

4 Mean-Field Stochastic Differential Equations 

We consider measurable functions b : fi x [0, T] x R n x R n -» R n and a : Q x [0, T] x M n xK n ^ 
jgmxd w j 1 j c ] 1 are supposed to satisfy the following conditions: 

(i) &(•, 0, 0) and cr(-, 0, 0) are F — progressively measurable continuous processes and there exists 
some constant C > such that 

\bit,x',x)\ + \ait,x',x)\ < C(l + \x\),a.s., for all < t < T, x,x' £R n ; 
(ii) b and a are Lipschitz in x, x' , i.e., there is some constant C > such that 

|6(t, x[,xi) - bit,x' 2 ,x 2 )\ + \ait,x[,xi) - ait,x' 2 ,x 2 )\ < Ci\x[ - x' 2 \ + \xi - x 2 \), a.s., 

for all < t < T, xi, x[, x 2 , x 2 G R n . 

(H4.1) 
We now study the following SDE parameterized by the initial condition (£, Q G [0, T] x 
L 2 in,F t ,P;R n ): 

j dX 1 / = E , [b(s,(x¥y,x¥)]ds + E , [<T( S ,(X t /) , ,x¥)]dB a , se[t,T], 

xl< = C- (41) 



Recall that 



E'[bis,iX t /y,X t /)]iu;)= / biJ,u,s,x¥iu>),x¥iu))Pidw'), ,e!]. 
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Theorem 4.1. Under the assumption (H4-1) SDE (4-1) has a unique strong solution. 

The proof is made in two steps like that of Theorem 3.1 and it uses standard arguments for forward 
SDEs. Since the proof is straight-forward we prefer to omit it. 

Remark 4.1. From standard arguments we also can get that, for any p > 2, there exists C p G R 
such that, for all t G [0,T] and (,(' G LP(Q, ^i,P;R n ), 

E[ sup \X 1 / - X*/' \P\F t ] < C P |C-CT, a.s., 

t<s<T 



E[ sup |X*' C |TO < C p (l + |Cn, a.s., ( 4. 2) 



t<s<T 

e[ sup |x*' c -cn^] < c p (i + icn<5i, 

t<s<t+<5 

P-a.s., for all 5 > iwii/i £ + <5 < T. 

These in the classical case well-known standard estimates can be consulted, for instance, in 
Ikeda, Watanabe J8$, pp. 166-168 and also in Karatzas, Shreve f^, pp. 289-290. We also emphasize 
that the constant C p in (4-2) only depends on the Lipschitz and the growth constants of b and a. 

5 Decoupled Mean-Field Forward-Backward SDE and Related DPP 

In this section we study a decoupled Mean-Field forward-backward SDE and its relation with 
PDEs. Given continuous functions b : U x [0, T] x R n x R n -> R n and a : Cl x [0, T] x R n x R n -> R nxd 
which are supposed to satisfy the conditions (i) and (ii) of (H4.1) and an arbitrary xo G R n , we 
consider the following SDE parameterized by the initial condition (£, C) G [0, T] x L 2 (£l, J-t, P;M. n ): 

f dX f / = E'[b(s, (X^Y^t^ds + E'[a(s, (#»)', X«)]dB a , s G [t,T], 

{ ^ = c . (5 ' 1} 

Under the assumption (H4.1), SDE (5.1) has a unique strong solution. Indeed, from Theorem 4.1 we 
first deduce the existence and uniqueness of the process X°' x ° G <S|(0, T; R n ) to the Mean-Field SDE 
(5.1). Once knowing X 0,x °, SDE (5.1) becomes a classical equation with the coefficients b(u>, s, x) = 
E'[b(uj', tu, s,X s ' x °(u)'),x)] and a(u>,s,x) = E'[a(uj',uj, s,X s ' x °(u'),x)]. Combining estimate (4.2) 
for (t, (") = (0, xo) with standard arguments for SDEs we obtain (4.2) also for equation (5.1). 

Let now be given two real-valued functions f(t,x',x,y',y,z) and <J>(x',x) which shall satisfy 
the following conditions: 

(i) $ : il x R n x R n — ► R is an Tt <8> B(M. n ) -measurable random variable and 

/ : n x [0, T] x R n x R n x R x R x R d -» R is a measurable process such that 
f{-,x',x,y',y,z) is F-adapted, for all (x',x,y f ,y,z) G R n x R n x R x R x R d ; 
(ii) There exists a constant C > such that 

\f(t,x' 1 ,x 1 ,y' 1 ,y 1 ,z 1 ) - f(t,x' 2 ,X2,y' 2 ,y2,Z2)\ + [$(a/ 1? xi) -$(x' 2 ,x 2 ) 

< Cdxi - x' 2 | + |xi - x 2 | + |yi - y' 2 \ + I2/1 - 2/2 1 + \z\ - z 2 \), a.s., 
for all < t < T, xi, x[, x 2 , x' 2 G R n , j/i, y[, y 2 , y' 2 G R and zi, z 2 G R rf ; 

(iii) / and $ satisfy a linear growth condition, i.e., there exists some C > 
such that, a.s., for all x', x G R n , 

|/(t,x',x,0,0,0)| + |$(x',x)| < C(l + \x\ + \x'\); (H5.1) 

(iv) f(Q,t,x',x,y',y,z) is continuous in £ for all (x',x,y',y,z), P(d(D)-a.s.; 
(v) /(£, x', x, y', y, z) is nondecreasing with respect to y' . 
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We consider the following BSDE: 

[ -dY* = E f \f(a,(x2*°)\xP,(X?™y,Y*'<,zF)]d8-ZpdB a ,8£[t,n ( , 9) 

| Yfi = E'[$((X% Xo y,X^)}. ( } 

We first consider the equation (5.2) for (t, £) = (0, xo): We know from Theorem 3.1 that there exists 
a unique solution (Y°> x °, Z°' x °) G S|(0, T; R) x W|(0, T; R d ) to the Mean-Field BSDE (5.2). Once we 
have (Y°' x °, Z 0,x °), equation (5.2) becomes a classical BSDE whose coefficients f(u),s,X 8 ,y,z) = 
E'[f(.,u, s, {Xs' X0 Y, X 1 /, (Y s °' xo y, y, z)] satisfies the assumptions (Al) and (A2), and $(w, X l ^{u)) = 
P'[<3?(., uj, (X T ' x °y, Xj, )] G L 2 (Q, J-j>, P). Thus, from Lemma 2.1 we know that there exists a unique 
solution (Y^,Z<^) G S|(0,T;R) x HJ(0,T;R d ) to equation (5.2). 

By combining classical BSDE estimates (see, e.g., Proposition 4.1 in Peng [J3]; or Proposition 
4.1 in El Karoui, Peng and Quenez [7]) with the techniques presented above we see that there 
exists a constant C such that, for any t G [0, T] and £, £' G L 2 (J7, ^j, P; R 71 ), we have the following 
estimates: 

(i)E[ sup |y'' c | 2 + jf |z|' c | 2 ds|.F t ] < C(l + |C| 2 ), a.s.; 

t<s<T 

(u)E[ sup lY^ - y s *' C '| 2 + f. T |Z*' C - Z*' C '| 2 ^|^] < C\( - C'| 2 , a.s. 
(iii) |y/' C |<C(l + |C|), a.s.; 

(iv) |y^-y«'|<qc-aa. s . (5 ' 3) 



In particular, 



Here the constant C > depends only on the Lipschitz and the growth constants of b, a, f and <J>. 
Let us now introduce the random field: 

u(t,x) = Y^ x \ s=t , (t,x) G [0,71 x R n , (5.4) 

where Y l > x is the solution of BSDE (5.2) with x G R n at the place of ( G L 2 (n, T u P;R n ). 

As a consequence of (5.3) we immediately have that, for all t G [0, T], P-a.s., 

(i) \u(t,x) — u(t,y)\ < C\x — y\, for all x,y G R n ; , , 

(ii) \u{t,x)\ < C(l + |x|), for all x G R n . ^ " ' 

Remark 5.1. In the general situation u is an adapted random function, that is, for any x G 
R n , u(-,x) is an F— adapted real-valued process. Indeed, recall that b,o~ and f all are F-adapted 
random functions while $ is Ft -measurable. However, if the functions b, a, f and <J> are determin- 
istic it is well known that also u is a deterministic function of (t,x) (see, e.g., Proposition 2.4 in 
Peng \T^). 

From now on, let us suppose that 

(vi) The coefficients b, a, f and <£ are deterministic, i.e., independent of . . 

(w',u;) efi = fi x n. ( ' ' 

The function u and the random field Y^, (t, () G [0,T] x L 2 (ft,.P t ,P;M n ), are related by 
the following theorem. 
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Theorem 5.1. Under the assumptions (H4-1) and (H5.1), for any t £ [0, T] andC, € L 2 (Q,JF i; P;R n ) 
we /iaue 

n(t,C) = ^/' C , P-o-s.. (5.6) 

Indeed, once X ^ and Y 0,x ° are determined, the coefficients of SDE (5.1) and BSDE (5.2) 
are well-determined, deterministic, depend only on (t,x) and (t,x,y,z), respectively, and satisfy 
the standard growth and Lipschitz conditions. Hence, Theorem 5.1 is a consequence of the corre- 
sponding result in Peng [TJ] (Theorem 4.7) (or, see also Theorem 6.1 in Buckdahn and Li [3]). 

We now discuss (the generalized) DPP for our FBSDE (5.1), (5.2). For this end we have to 
define the family of (backward) semigroups associated with BSDE (5.2). This notion of stochastic 
backward semigroup was first introduced by Peng [J3] and originally applied to study the DPP for 
stochastic control problems. Our approach extends Peng's ideas to the framework of Mean-Field 
FBSDE. However, we change the definition of the stochastic backward semigroup to simplify the 
proof of the existence of a viscosity solution of the associated PDE. 

Given the initial data (t,x), a positive number 5 <T — t and a real- valued random variable 
r, E L 2 (0,J^ +<5 ,P;M), we put 

G^ t+ M--=Ys' x , se[t,t + 5], (5.7) 

where the couple (Y s ' x , Z s ' x ) t < s < t+ s is the solution of the following BSDE with the time horizon 
t + 8: 

j -(^* = ^[/( a) (# o )^^^(y^ o )^lf^z*•*)]d S -z2• a! dB S) se%t+s], ._ a . 

1 Y t+5 = m 

here X t,x is the solution of SDE (5.1) and Y t,x is the solution of BSDE (5.2). Then, obviously, for 
the solution (Y*»*, Z^ x ) of BSDE (5.2) we have 

Gg.[*(JT^)] = Gg^pjg], < t < t + 6 < T. (5.9) 

Let us point out that in difference to Peng's definition of the backward semigroup the driver 
of our BSDE depends on the processes Y 0,x ° and Y t,x given by BSDE (5.2). The choice to let the 
driver / depend on these processes and not on Y 0,x ° , Y l ' x will simplify the proof of the existence 
theorem for the associated nonlocal PDEs in Section 6. 

Moreover, we have the following DPP 



u(t,x) = Yf* = G^MXp)} = G% +s [Y t f s ] 



t,t+5^t+5\ (510) 



whose simple form explains by the fact that our stochastic evolution system doesn't depend on a 
control. Here, for the latter relation we have used, that due to the uniqueness of the solution of 



rt.X 



SDE (5.1) and of BSDE (5.2), Y^ 5 = Y^' Xt ' +s , P-a.s., so that then from Theorem 5.1 it follows 
that Y t f s = u(t + 5, X t f s ), P-a.s. 

Remark 5.2. If f is independent of(y,z) it holds that 

-t+8 



G t £ +s [r ] }=E[r ] + [ E'[f(r, (*?-*»)', *£*, (Y r ' x °)')]dr\T s l se[t,t + 8], P- 

J s 



a.s. 
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In (5.5) we have already seen that the value function u(t, x) is Lipschitz continuous in x, 
uniformly in t. Relation (5.10) now allows also to study the continuity property of u(t,x) in t. 

Theorem 5.2. Let us suppose that the assumptions (HJ^.l), (H5.1) and (H5.2) hold. Then the 
value function u(t, x) is -^ — Holder continuous in t, locally uniformly with respect to x: There exists 
a constant C such that, for every x € M. n , t, t' £ [0, T], 

\u(t,x) -u(t',x)\ < C(l + \x\)\t-t'\K 

Proof. Let (t, x) € [0, T] x W l and 5 > be arbitrarily given such that < 5 < T — t. Our 
objective is to prove the following inequality by using DPP: 

-C(l + \x\)6v <u(t,x) -u(t + S,x) < C(l + |x|)<$3. (5.11) 

From it we obtain immediately that u is ^—Holder continuous in t. We will only check the second 
inequality in (5.11), the first one can be shown in a similar way. To this end we note that due to 
(5.10), 

u(t,x)-u(t + 5,x)=l}+l], (5.12) 

where 

I] := Gf t+s [u{t + 6,X% 5 )]-Gf t+8 [u{t + ^)l 

l l ■= G^ t+s [u{t + 5,x)]-u{t + 5,x). 

From Lemmata 2.3 and 3.1 and the estimate (5.5) we obtain that, 

\I}\ < [CE(\u(t + 6,X% s )-u(t + 5,x)\ 2 \f t )} 1 * 
< [CE(\X^ 5 -x\^ t )]K 

and since E[\X t ' x s — x| 2 |jT t ] < (7(1 + \x\ 2 )5 we deduce that \lj\ < C(l + \x\)5z. From the definition 
of (7 t '^ +( j[-] (see (5.7)) we know that the second term L 2 can be written as 

I 2 = E[ U (t + S,x) + J t t+S E'[f( S ,(X^ Xo )\xi' x ,(Y s °' Xo yX'^Zt' x )}ds 
- J t Z s ,x dB s \F t ] - u(t + S, x) 

= E[f t t+5 E'lf(s,(X°> Xo y,Xi> x ,(Y s °^y,Y s ^ \Z^)]ds\T t ]. 

Then, with the help of the Schwartz inequality, and the estimates (4.2), (5.3)-(i) for the BSDEs 
(5.2) and (5.8) (with n = u(t + 5, x)) and (5.5) we have 

|/|| <ShE[f^\E , [f( S ,(X^)'X^,(Y°^y,Y^Jl' x )]\ 2 ds\^ t ]h 

<(y3£[/^(|tf[/( s ,(x, ^ 

< CS^E[f^\\l + \xt> x \ + \Yt x \ + \zt' x \) 2 ds\T t ]l 
<C{l + \x\)5^. 

Hence, from (5.12), we get the second inequality of (5.11) 

u(t,x) -u(t + 5,x) < (7(1 + \x\)8*. 
The proof is complete. D 
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6 Viscosity Solution Of PDE: Existence Theorem 

In this section we consider the following PDE 

gu(t, x) + Au(t, x) + E[f(t, Xt' xo ,x, u(t, X^ x °),u(t, x), Du(t, x).E[a(t, X^ xo , x)])] = 0, 

(t,x) E [0,T) xl", 
u(T,x) = E[<S>(X% x °,x)], x e K n , 



(6.1) 



with 



Au(t, x) : = hr{E[a{t, X° t ' x ° , x)]E[a(t, X^ x ° , x)] T D 2 u(t, x)) + Du(t, x).E[b(t, X^ x ° , x)]. 

Here the functions b,a,f and $ are supposed to satisfy (H4.1), (H5.1) and (H5.2), respectively, 
and X°' xo is the solution of the Mean-Field SDE (5.1). 

We attract the reader's attention to the fact that, since 

E[f(t, X^ , x, u(t, X?' X0 ),u(t, x),Du(t, x).E[a(t, X^ x ° , x)})} 
= Jr" /(*> x ' ■> x > u (*> x ')i u (^> x),Du(t, x).E[a(t, X t ,x ° , x)])P x o,x (dx'), 

the above equation is in fact a nonlocal PDE. 

In this section we want to prove that the value function u(t, x) introduced by (5.4) is the 
viscosity solution of equation (6.1). For this we extend Peng's BSDE approach [T3] developed 
in the framework of stochastic control theory to that of the Mean-Field FBSDE. The difficulties 
related with this extension come from the fact that now, contrarily to the framework of stochastic 
control theory studied by Peng, we have to do with nonlocal PDEs. Moreover, in difference to [3] 
the nonlocal term is not generated by a diffusion process with jumps. This fact is the source of 
difficulties mainly in the proof of the uniqueness of the viscosity solution (given in the next section) 
which are different to those in [3] . Let us first recall the definition of a viscosity solution of equation 
(6.1). The reader more interested in viscosity solutions is referred to Crandall, Ishii and Lions |6j. 

Definition 6.1. A real-valued continuous function u £ C p ([0,T] x W 1 ) is called 
(i) a viscosity subsolution of equation (6.1) if, firstly, u(T,x) < E[$(X T ' X0 ,x)], for all x £ M. n , and 
if, secondly, for all functions if £ Cf b ([0, T] x M. n ) and (t,x) £ [0, T) x W 1 such that u — if attains 
its local maximum at (t,x), 

&(p(t, x) + D<p(t, x).E[b(t, X° t ' x ° , x)\ + \tr(E[o-{t, X?' x ° , x)]E[a(t, X?' x ° , x)\ T D 2 <p{t, x)) 
+E[f(t, X°' xo , x, u(t, X?> xo ),u(t, x),D<p{t, x).E[a(t, X^ x ° , x)\)\ > 0; 

(ii) a viscosity supersolution of equation (6.1) if, firstly, u(T,x) > E[$>(X T ' X ° ,x)], for all x € W 1 , 
and if, secondly, for all functions tp £ Cf b ([0,T] x M. n ) and (t,x) £ [0, T) x M. n such that u — ip 
attains its local minimum at (t,x), 

%<p(t, x) + D<p(t, x).E[b(t, X?> x ° , x)] + ltr(E[a(t, X?' x ° , x)]E[a(t, X^ , x)} T D 2 <p(t, x)) 
+E[f(t,X t , ' Xo ,x,u(t,X t ) > Xo ),u(t,x),D<p(t,x).E[o-(t,X i , ' Xo ,x)])]<0; 

(hi) a viscosity solution of equation (6.1) if it is both a viscosity sub- and a supersolution of equation 
(6.1). 
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Remark 6.1. (i) C p ([0,T] xR n ) = {u € C([0, T]xR n ) : There exists some constant p > such that 

\u(t,x)\ i -i 

su P(t,x)e[0,T]xi?" 1+|i|p < +°°J"- 

(ii) Usually the definition of a viscosity solution is given for test fields ip of class C 1,2 . However, 
it can be shown that it is sufficient to work with test functions from Cf b ([0, T] x R n ). The space 
Cf b ([0,T] x M. n ) denotes the set of the real-valued functions that are continuously differentiable up 
to the third order and whose derivatives of order from 1 to 3 are bounded. 

We now can state the main statement of this section. 

Theorem 6.1. Under the assumptions (H^.l), (H5.1) and (H5.2) the function u(t,x) defined by 
(5.4) is a viscosity solution of equation (6.1). 

The proof of the theorem uses the BSDE method of Peng [T3]. However, it is simplified by 
the specific choice of our stochastic backward semigroup. For the proof of this theorem we need 
four auxiliary lemmata. To abbreviate notations we put, for some arbitrarily chosen but fixed 
V eCf >b ([0,T]xM n ), 

F(s, x, z) = §- s y{s, x) + \tr(da T {s, x)D 2 (p) + Dtp.b(s, x) , . 

+f u (s,x,z + Dtp(s,x).a(s,x)), [ ' ' 



where 



a(s,x) = E[o-(s,X°' X0 ,x)},b(s,x) = E[b(s,X°> X0 ,x)}; 
f u (s,x,z) = E[f(s,X°' X0 ,x,u(s,X°> X0 ),u(s,x),z)}, 



(s, x, z) e [0, T]xt n x M. d , and we consider the following BSDE defined on the interval [£,< + £] (0 < 
5<T-t): 

f -dY} = F{s, Xt' x , Z])ds - Z]dB a , 

1 

't+5 
where the process X t,x has been introduced by equation (5.1). 



Y i _ (6-3) 



Remark 6.2. It's not hard to check that F(s, X s ' x , z) satisfies (Al) and (A2). Thus, due to Lemma 
2.1 equation (6.3) has a unique solution. 

We can characterize the solution process Y 1 as follows: 
Lemma 6.1. For every s E [t, t + 5], we have the following relationship: 

Y s = G% +S W + 5, X% 5 )] - <p(s, Xl>% P-a.s. (6.4) 

Proof. We recall that Gj x +S [ip(t + 5, X t ' x s )] is defined with the help of the solution of the BSDE 

f -dY s = E'[f(s, (X^ X0 Y, Xt' x , (Y s °> xo )', Y*' x , Z s )]ds - Z s dB s , se[t,t + 5], 
\ Y t+5 =<p(t + 5,X% s ), 

by the following formula: 

G% 5 W + S, X% 6 )] =Y S , se [t, t + 5] (6.5) 
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(see (5.7)). We also recall that E , [f(s,(X%' aso y,X t a ' x ,(Y?' xo y,Yt' x ,Z a )] = f u (s,Xi' x ,Z S ). There- 
fore, we only need to prove that Y s — (p(s, X/ x ) = Y^. This result can be obtained easily by applying 
Ito's formula to (p(s,X s ,x ). Indeed, we get that the stochastic differentials of Y s — (p(s,X/ x ) and 
l^, 1 coincide, while at the terminal time t + S, Y+s — f{t + 5,X t ' x s ) = = Y^ +s . So the proof is 
complete. □ 

We now introduce the deterministic function 

i-t+6 



/i-t-o 
F(r,x,0)dr, s€[t,t + S\. 



Obviously, the couple (Y 2 ,Z 2 ) = (Y 2 ,0) is the unique solution of the following (deterministic) 
BSDE in which the driving process X t,x is replaced by its deterministic initial value x: 



-dY s 2 = F(s,x, Z 2 )ds - Z 2 dB s , 
Y 2 +s =0, s£[t,t + 5}. 



2 n _ r, , , n ( 6 - 6 ) 



The following lemma will allow us to neglect the difference \Y^- — Y 2 \ for sufficiently small 5 > 0. 

Lemma 6.2. We have 

\Y?-Y?\ <C6%, P-a.s. (6.7) 

Proof. We recall that from (4.2) with £ = x it follows that there is some constant C p depending 
on p, x, but not on 8 > 0, such that 

E[ sup \Xl> x - x\ p \f t ] < CpS* . (6.8) 

t<s<t+S 

We now apply Lemma 2.3 combined with (6.8) to the equations (6.3) and (6.6). For this we set in 
Lemma 2.3: 

Ci =6 = 0, g(s,y,z) =g(s,z) = F (s , X f s ' x , z) , 

¥>i(s) = 0, M^ = F(s,x,Z 2 ) - F{s,X l s ' x ,Z 2 s ). 

Obviously, the function g is Lipschitz with respect to z, and |</?2(s)| < C(l H- |cc| 2 )(|-X" s ,a; — x| + |X s ' x — 
x\ 3 ), for s G [t,t + S],(t,x) € [0,T) xR". Thus, with the notation p (r) = (1 + |x| 2 )(r + r 3 ), r > 0, 

we have 

E[jt +S \Zl - Z 2 \ 2 ds\f t ] 

< CE[J t t+S pl{\xl' x -x\)ds\F t ] 

< C5E[ sup p 2 (\xi' x - x\)\T t ] 



t<s<t+S 

< C5 2 . 



Therefore, 



\Y?-Y?\ = \E[0?-Y?)\r t ]\ 
= \E[J t t+S (F( S ,Xi' x ,Zl) - F(s,x,Z 2 ))ds\F t }\ 

< CE[J t t+S [p (\X^ x -x\) + \Z] - Z 2 \]ds\T t ] 

< CE[f t t+5 p (\xi> x - x\)ds\T t ] + CS^E[J t t+6 \Z l s - Z 2 \ 2 ds\^ 



< C8 3 2. 
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Thus, the proof is complete. □ 

Now we are able to give the proof of Theorem 6.1: 
Proof. Obviously, u(T, x) = E[^(Xj! x ° , x)], x G M. n . Let us show that u is a viscosity supersolution 
(respectively, subsolution). For this we suppose that <p G Cf b ([0,T] x W 1 ) and (t,x) G [0, T) x R n 
are such that u — ip attains its minimum (respectively, maximum) at (t,x). Notice that we can 
replace the condition of a local minimum (respectively, maximum) by that of a global one in the 
definition of the viscosity supersolution (respectively, subsolution) since u is continuous and of at 
most linear growth. Moreover, without loss of generality we may also suppose that </?(£, x) = u(t, x). 
Then, due to the DPP (see (5.10)), 

<p(t,x) = u(t, x) = G^ t+S [u{t + 8,X%s)], 0<5<T-t, 

and from u > (p (respectively, u < tp) and the monotonicity property of G t ' x , ?[■] (see Lemma 2.2 
and Theorem 3.2) we obtain 

G t 't +S [(p(t + 5,X t '^ s )] — (p(t,x) < (respectively, > 0), P-a.s. 

Thus, from Lemma 6.1, 

Y t < (respectively, > 0), P-a.s., 

and furthermore, from Lemma 6.2 we have 

ft+S 3 3 

/ F(s,x,0)ds = Y t 2 < C5* (respectively, > -C8^), P-a.s. 

It then follows that 

F(t, x, 0) < (respectively, > 0) 

and from the definition of F we see that u is a viscosity supersolution (respectively, subsolution) 
of equation (6.1). Finally, we prove that u is a viscosity solution of equation (6.1). □ 

7 Viscosity Solution of PDE: Uniqueness Theorem 

The objective of this section is to study the uniqueness of the viscosity solution of PDE (6.1), 

§- t u(t, x) + Au(t, x) + E[f(t, Xt' xo ,x, u(t, Xt' xo ),u(t, x),Du(t, x).E[a{t, X^ xo , x)])] = 0, 

(t,x) e[0,T) xR n , 
u(T,x) = E[<S>(X^ x °,x)], x G W 1 , 

(7.1) 
with 

Au(t, x) := Ur(E[a(t, X°' x ° , x)]E[a(t, X^ x ° , x)] T D 2 u{t, x)) + Du(t, x).E[b(t, X^ x ° , x)]. 

Here the functions b,a,f and $ are supposed to satisfy (H4.1), (H5.1) and (H5.2), respectively. 

We will prove the uniqueness for equation (7.1) in the space C p ([0, T] x W n ) of continuous 
functions with at most polynomial growth. In an earlier work Barles, Buckdahn and Pardoux [3] 
introduced the space of continuous functions 
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e = {^ G C([0,T] xi n ):]4>0 such that lim^^ \ip(t,x)\ exp{-A[log((|:r| 2 + l)?)] 2 } = 
0, uniformly in t G [0, T]}. 

Its growth condition is slightly weaker than the assumption of polynomial growth but more restric- 
tive than that of exponential growth. They proved in the uniqueness of the viscosity solution 
of an integro-partial differential equation associated with a decoupled FBSDE with jumps. It was 
shown in [3] that this kind of growth condition is optimal for the uniqueness. However, as the 
following example shows, we cannot hope to have this class for the uniqueness also for our type 
of PDE. 

Example 7.1 Let n = d = 1, o~(s,x',x) = ax, b(s,x',x) = ^-x and xq = 1. Then X s ' x ° = 
exp{o\Bg}, s G [0, T], and for f(s,x',x,y',y,z) = y' , PDE (7.1) takes the form 

f £ t u(t, x) + %x&u(t, x) + ^-x 2 £ju{t, x) + E[u{t, exp{aB t })] = 0, {t, x) G [0, T) x R, 
y u(T, x) = E[$(ex.-p{aB T }, x)], x G R. 

But, for the growth function 

X(t,x) = exp{i[log((|x| 2 + l)i)] 2 }, (t,x) G [0,T) x R, 

which comes from the definition of 0, we have 

E[x(t,X^ )]= E[exp{A[log((\exp{aB t }\ 2 + 1)^}} 
> E[exp{Aa 2 B?}} 

= +oo, iftG[^,T]. 

As the example shows, for a function u G 0, the coefficient E[f(t,X t ,X0 ,x,u(t,X t ' x °),y,z)] may 
be not well defined. This is the reason why we restrict the study of the uniqueness to the smaller 
class C p {[0,T] xl"). 

Theorem 7.1. We assume that (H4.I), (H5.1) and (H5.2) hold. Let u\ (resp., u 2 ) G C p ([0,T] x 
R n ) be a viscosity subsolution (resp., supersolution) of equation (7.1). Then we have 

ui{t,x) < u 2 (t,x), for all (t,x) G [0,T] x R n . 

The proof of the theorem will be prepared by the following auxiliary lemmata. 

Lemma 7.1. Let K be a Lipschitz constant of f{t,x' , ., ., ., .), uniformly in (t,x') from (H5.1)-{\\), 
and let v > K. Then, if u G C p ([0, T] x R n ) is a viscosity subsolution (respectively, supersolu- 
tion) of PDE (7.1) the function u(t,x) = u{t,x)e ut , (t,x) G [0,T] x M n , is a viscosity subsolution 
(respectively, supersolution) of the following PDE: 

gu(t, x) + Du(t, x).E[b(t, X^ xo , x)} + \tr{E[a{t, X^ x ° , x)]E[a(t, X°' x ° , x)] T D 2 u(t, x)) 
+E[f{t,X^ Xo ,x,u{t,X^ X0 ),u(t,x),Du(t,x).E[a(t,X^ Xo ,x)})} =0, (t,x) G [0, T) x R n ; 
u(T,x) = E[$>{X^ x °,x)]e vT , x G R", 

(7.2) 
where f(t, x', x, y', y, z) = e ut f(t, x', x, e~ ut y', e~ ut y, e~ ut z)-vy, (t, x' , x, y',y, z) G [0, T] x R n x R n x 
R x R x R rf , conserves the properties (H5.1) and (H5.2) of f and is, moreover, strictly decreasing 
in y: 

f{t,x',x,y',yi,z) - f(t,x',x,y',y 2 ,z) < -{y - K)(yx -y 2 ), 

for all (t,x',x,y', z)and all y±,y 2 G R with y\ > y 2 . 
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The proof of this well-known transformation is straight-forward and is hence omitted. 

Lemma 7.2. Let u\ G C p ([0,T] x R n ) be a viscosity subsolution and u^ G C p ([0, T] x R n ) be a 
viscosity super solution of equation (7.2). Then the function to := u\ — u<i G C p ([0,T] x W l ) is a 
viscosity subsolution of the equation 

-vu(t,x) + §- t uj(t,x) + \tr(aa T (t,x)D 2 uj) +Du>.b(t,x) + K\to(t,x)\ + 

+ KE[(Lu{t,Xt> Xo )) + }+K\Duj.a{t,x)\=0, (t,x) G [0,T) x R n , (7.3) 

u(T,x) = o, ier. 

(See (6.2) for the notations b and a.) 

The proof of this lemma follows from that of Lemma 3.7 in [3], it turns out to be even simpler 
because contrary to Lemma 3.7 in [3] we don't have any integral part generated by a diffusion 
process with jumps here in the equations (5.1) and (5.2). 

Lemma 7.3. For any A > 1, there exists C\ > such that the function 

X (t,x) = Ae c ^ T - t ^(x), 

with 

tj){x) = (|x[ 2 + l)i, x GM n , p> 1, 

satisfies 



§- tX (t, x) + \tr{da T {t, x)D 2 X ) + D X ~b{t, x) + K X (t, x) + 
K\Dx{t,x).a{t,x)\+KE[x{t,Xt' X0 )} < in [0,T] x K n 



(7.4) 



Proof. By direct computation we first deduce the following estimates for the first and second 
derivatives of i/j: 

\d^x)\ = P -^-\x\, |z?V(*)l<P 2 nrr7> ^ Rn ' 

\xr + 1 \xr + 1 



and 



d - X (t,x) = -C lX (t,x), E[i,(X^ X0 )]<C p i;(x), xG 



dt' 

The latter estimate is a direct consequence of (4.2). Taking into account that the coefficients b, a 
are of most linear growth, these estimates imply for all t G [0, T], 

| X (t, x) + ±tr(aa T (t, x)D 2 X ) + D X &(t, x) + K X (t, x)+ 

K\D X (t, x).a(t, x) | + KE[ X (t, X?> x °)\ 
< - X (t,x){d - p 2 C - K - KC P } 
= -x(t, x) < 0, if C\ := p 2 C + K + KC P + 1 large enough. 

□ 

Now we can prove the uniqueness theorem. 

Proof of Theorem 7.1. From Lemma 7.1 we only need to prove that: If u\ (respectively, U2) 

G C p ([0,T] x M n ) is a viscosity subsolution (respectively, supersolution) of equation (7.2) then we 

have 

ui(t,x)<u 2 (t,x), for all (t,x) G [0, T] x R n . 
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Let us put (jj := u\ — v,2- Since to G C p ([0,T] x M n ), there exist some constants C > and q > 1 
such that \to(t,x)\ < C(l + [x|) 9 , (i,x) € [0, T] x M n . In the definition of 

X (t,x) = Ae Cl{T - t] {\x\ 2 + 1)5 

we now choose p > q. Then, for any a > 0, u(t, x) — ax(t, x) is bounded from above in [0, T] x R n 

and 

M := max (u - ax)(t, x)e~ K{ - T ' t) 
[0,T]xR™ 

is achieved at some point (ioi^o) € [0, T] x R n (depending on a). We now have to distinguish 

between two cases: 

For the first case we suppose that u>(to,xo) < 0, for all a > 0. 

Then, obviously M < and u±(t,x) — U2(t,x) < ax(t,x) in [0,T] x M n . Consequently, letting a 

tend to zero we obtain 

ui(t,x) < u 2 (t,x), for all (t,x) € [0, T] x M n . 

For the second case we assume that there exists some a > such that u>(to, xq) > 0. 

We notice that oo(t, x) - ax(t, x) < (co(t ,x ) - ax(t , x ))e- K( - t - t °^ in [0, T] x W 1 . Then, putting 

<p(t, x) = a X (t, x) + (u- a X ){to,xo)e- K ^- to \ (t, x) G [0, T] x R n , 

we get to — ip < = (lo — (p)(to,xo) in [0, T] x R n . Consequently, since uj is a viscosity subsolution 
of (7.3) from Lemma 7.2, 



-vu(t , x ) + $t<p{t , x ) + 2tr(a<J T (to,x )D 2 (p(to, x )) + £>¥>(*o, a; ).6(t , ^o)+ 

0,3 

to 



K\<p(t Q ,x )\ + K\D<p(t Q ,x Q ).*(t Q ,xo)\ + KE[(^to,X?f°))+} > 0. 



Moreover, due to our assumption that u(to,Xo) > and since to(to,xo) = ip(to,xo) we can replace 
K\(p(to,xo)\ by Kip(to,xo) in the above formula. Then, from the definition of (p and Lemma 7.3, 

< -vLu(t , x ) + m^tto, xo) + \tr(aa T (tQ, x )D 2 (p(t , x )) + D(p(t , x ).b(t , x )+ 

K\cp(t ,x )\ + K\Dip(t ,xo).a(to,x )\ + KE[(<p(t ,X% X0 ))+] 
= a^x(*o, zo) + "^X^o, »o) + a±tr(o-<5- T (i , a; )Z? 2 x(*o, ^o)) + a£>x(*o, »o)-K*o, £o) + 

aif|Ux(*o,a;o)-o-(*o,a;o)| + KE[(ax(to,X^ x °) + (w - ax)(£o,^o)) + ] - w(t ,zo) 

\ 



< a{^{t ,x ) +Kx{t ,x ) + \tr{aa T \t , x )D 2 x(t , x )) + Dx(to,x ).b(to,x ) 

+K\Dx(to,xo).Hto,xo)\ + KE[ x (to,X?f )}} < 0, 



since v > K, which is a contradiction. Thus, the proof is complete. □ 

Remark 7.1. Obviously, since the value function u(t,x) is of at most linear growth it belongs to 
C p ([0,T] x R n ), and so u(t,x) is the unique viscosity solution in C p ([0,T] x M. n ) of equation (6.1). 
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